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Using the light-front quark model, we calculate the 
transition form factors, decay rates, and longitudinal lep- 
ton polarization asymmetries for the exclusive rare B 
K£'^£~{£ = e,ii,T) decays within the standard model. Eval- 
uating the timelike form factors, we use the analytic continu- 
ation method in q"*" = frame to obtain the form factors F+ 
and Ft, which are free from zero- mode. The form factor F- 
which is not free from zero-mode in = frame and con- 
taminated by the higher(or nonvalence) Fock states in g"*" 7^ 
frame is obtained from an effective treatment for handling the 
nonvalence contribution based on the Bethe-Salpeter formal- 
ism. The covariance(i.e. frame-independence) of our model 
calculation is discussed. We obtain the branching ratios for 
BR(B ^ K£+r) as 4.96 x 10-'^\Vts/Vcb\'^ for £ = e,fj. and 
1.27 X 10-''\Vts/V,b\^ for £ = t. 



I. INTRODUCTION 

The upcoming and currently operating B factories 
BaBar at SLAC, Belle at KEK, LHCB at CERN and 
B-TeV at Fermilab as well as the planned r-Charm fac- 
tory CLEO at Cornell make the precision test of standard 
model(SM) and beyond SM ever more promising Es- 
pecially, a stringent test on the unitarity of Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrix in SM will 
be made by these facilities. Accurate analyses of exclu- 
sive semileptonic B-decays as well as rare B-decays are 
thus strongly demanded for such precision tests. One 
of the physics programs at the B factories is the exclu- 
sive rare B decays induced by the flavor-changing neutral 
current(FCNC) transition. Since in the standard model 
they are forbidden at tree level and occur at the lowest 
order only through one-loop (Penguin) diagrams |§-D, 
the rare B decays are well suited to test the SM and 
search for physics beyond the SM. While the experimen- 
tal tests of exclusive decays are much easier than those of 
inclusive ones, the theoretical understanding of exlcusive 
decays is complicated mainly due to the nonperturbative 
hadronic form factors entered in the long distance non- 
perturbative contributions. The calculations of hadronic 
form factors for rare B decays have been investigated by 
various theoretical approaches, such as relativistic quark 
model [^-[1^1 , heavy quark theory ||ll[], three point QCD 
sum rules [|12[, light cone QCD sum rule [[l3| - [l6| , and 
chiral perturbation theory [|l7|JlS]. Perhaps, one of the 



most well-suited formulations for the analysis of exclu- 
sive processes involving hadrons may be provided in the 
framework of light-front quantization p9| . 

The aim of the present work is to calculate the hadronic 
form factors, decay rates and the longitudinal lepton po- 
larization asymmetries for B —> K£~^£~{£ = e,/i, and 
r) decays within the framework of the SM, using our 
light-front constituent quark model(LFCQM or simply 
LFQM)) based on the LF quantization. The lon- 

gitudinal lepton polarization, as another parity-violating 
observable, is an important asymmetry |2^ and could be 
measured by the above mentioned B factories. In partic- 
ular, the T channel would be more accessible experimen- 
tally than e- or /x-channels since the lepton polarization 
asymmetries in the SM are known to be proportional to 
the lepton mass. Although some recent works have 
studied the lepton polarizations using the general form 
of the effective Hamiltonian including all possible forms 
of interactions, we shall analyze them within the SM as 
many others did. 

Our LFQM |20-2^] used in the present analysis has 
several salient features compared to other LFQM [^|j|] 
analysis: (1) We have implemented the variational prin- 
ciple to the QCD motivated effective LF Hamiltonian 
to enable us to analyze the meson mass spectra as well 
as various wavefunction-related observables such as de- 
cay constants, electromagnetic form factors of mesons in 



spacelike (q < 0) region 



(2) We have performed 



the analytical continuation of the weak form factors from 
spacelike region to the entire (physical) timelike region to 
obtain the weak form factors for the exclusive semilep- 
tonic decays of pseudoscalar mesons [Q. (3) We have 
recently presented in |Q an effective treatment of han- 
dling the higher Fock state (or nonvalence) contribution 
to the weak form factor in (7+ > frames, based on the 
Bethe-Salpetcr(BS) formalism (see also p3||). 

The explicit demonstration of our analytic continua- 
tion method using the exactly solvable model of (3 -I- 
l)-dimensional scalar fleld theory model can be found 
in @. The Drell-Yan-West {q^=q°+q^=0) frame is use- 
ful because only valence contributions are needed as far 
as the -component of the current is used. Our an- 
alytic solution in the q~^=0 frame as a direct applica- 
tion to the timelike region differs from the method used 
in [^js] where the authors used a simple parametric for- 
mula extracted from the small q^ behavior of a form fac- 
tor. However, some of the form factors in timelike exclu- 
sive processes receive higher Fock state contributions(i.e. 
zero-mode in q'^ — frame or nonvalence contribution 



1 



in g+ 7^ frame) within the framework of LF quantiza- 
tion. Thus, it is necessary to include either zero-mode 
contribution(if working in = frame) or the nonva- 
lence contribution (if working in 7^ frame) to obtain 
such form factors. Specifically, in the present analysis 
of exclusive rare B — > Ki^i" decays, three independent 
hadronic form factors, i.e. F+{q^), F_(q^) from the V- 
yl(vector-axial vector) current, and FT^q^) from the ten- 
sor current, are needed. While the two form factors F+ 
and Ft can be obtained from only valence contribution in 
q+ = frame without encountering the zero-mode com- 
plication it is necessary to include the nonvalence 
contribution for the calculation of the form factor F_. 
Our effective method of calculating novalence con- 
tributions has been shown to be quite reliable by check- 
ing the covariance of the model. Thus, we utilize both 
the analytic method in g+ = frame to obtain Ft) 
and the effective method in (7+ > frame to obtain F_, 
respectively. 

The paper is organized as follows. In Sec. II, we discuss 
the standard model effective Hamiltonian for the exclu- 
sive rare B — > K£^£~ decays and reproduce the QCD 
Wilson coefficients necessary in our analysis. The formu- 
las of the hadronic form factors, differential decay rates, 
and the longitudinal lepton polarization asymmetries are 
also introduced in this section. In Sec. Ill, we calculate 
the weak form factors F+((?^), F_(g^) and -Pr('?^) using 
our LFQM. To obtain F+{q'^) and FT{q^), we use the 
g+ = frame (i.e. q^ = —q^ < 0) and then analyti- 
cally continue the results to the timelike q^ > region 
by changing q± to iq± in the form factors. The form fac- 
tor F-{q'^) is obtained from our effective method in 
purely longitudinal <?+ > frames (i.e. q^ = q~^q~^ > 0). 
In Sec. IV, our numerical results, i.e. the form fac- 
tors, decay rates, and the longitudinal lepton polarization 
asymmetries for B K£^£~ decays, are presented and 
compared with the experimental data as well as other 
theoretical results. Summary and discussion of our main 
results follow in Sec. V. In the Appendix A, we list the 
QCD Wilson coefficients necessary for the rare B ^ K 
transition. In the Appendix B, we show the derivation of 
the differential decay rate for B K£'^£^ in the case of 
nonzero lepton(mf ^ 0) mass. In Appendix C, we show 
the generic form of our analytic solutions for the weak 
form factors in timelike region. 



II. OVERVIEW OF EFFECTIVE HAMILTONIAN 
IN OPERATOR BASIS 

The rare b s£~^£~ decay process can be represented 
in terms of the Wilson coefficients of the effective Hamil- 
tonian obtained after integrating out the heavy top quark 
and the bosons [01, i.e. 



V2 



'ts / , 



a(A)Oz(/i), 



(1) 



where is the Fermi constant, Vij are the CKM ma- 
trix elements and Ci(/i) are the Wilson coefficients. It is 
known that the Wilson coefficients C3 — Cq of QCD pen- 
guin operators O3 — Oq are small enough to be neglected 
and also the operator Og (^ G"^ , strong interaction field 
strength tensor) does not contribute to 6 — > s£'^£~ tran- 
sition. Thus, the relevant basis operators Oi{jl) to the 
rare h s£~^£~ decay are 



O2 
O7 

O9 



167r2 
167r2 
167r- 



mb{saCr^yPi^ba)F^^'' , 
■{s^rPLbo.){h^-fd). 



(2) 



where -Pj,(/j) = (l=F75)/2 is the chiral projection operator 
and F^'^ is the electromagnetic interaction field strength 
tensor. The Lorentz and color indices are denoted as 
/i(and v) and a(and 8), respectively. The renormaliza- 
tion scale fi in Eq. (|l|) is usually chosen to be /2 ~ rrib 
in order to avoid large logarithms, ln(Mvy/mf,), in the 
matrix elements of the operators Oi. The Wilson coeffi- 
cients Ci{mh) determined by the renormalization group 
equations(RGE) from the perturbative values Ci(Mw) 
are given in the literature(see, for example [^jQ). 

Since the operators Oi and O2 contribute to 6 ^ s£^£^ 
through cc-loops which again couple to £^£r through vir- 
tual photon, they can be incorporated into an "effective" 
Oq . The resulting effective Hamiltonian in Eq. (|^) has the 
following structure (neglecting the strange quark mass) 



v*v 



2 sctm''? Pbp£i'^£ 



+ Cf{mb)sj^PLb£r£ + Cw{mb)s-i^PLb£nd 



(3) 

The effective Wilson coefficient Cf^{s=q^/m'l) is given 
by [pPg 



Gf (s)^Gf (j) + yLD(.s), 

= ^9 (1 + ^^(^)) + YsT>{s) + Fld(s), (4) 

where the function 1sd(s) is the one-loop matrix element 
of O9, yLD(s) describes the long distance contributions 
due to the charmonium vector J/i/'j "^'i ' ' ' resonances via 
B K{J/ip,ilj' , • • •) ^ K£'^£~, and lu{s) represents the 
one-gluon correction to the matrix element of O9. Their 
explicit forms are given in the literature P,^,p8|-^ and 
also in the Appendix A of this work. For the numerical 
values of the Wilson coefficients and relevant parame- 
ters in obtaining Eq. (|^), we use the results given by 
Refs. p9|,p0|l: m* = 175 GeV, mt = 4.8 GeV, = 1.4 
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GeV, as{Mw) = 0.12, a,{rrn,) = 0.22, Ci = -0.26 
C2 = 1.11, C3 = 0.01, C4 = -0.03, C5 = 0.008 
Cg = -0.03, C7 = -0.32, Cg = 4.26, and Cio = -4.62. 

In Fig. we plot the effective Wilson coefficient Cg 
as a function of s. As the real part of Cg^, the thick(thin) 
solid line represents the result with(without) LD contri- 
bution, i.e. Yie{Cf){Re{Cf )). The imaginary (dotted 
line) part of Cg^ is the result without LD contribution, 
lni(C'§^). In our numerical calculation of Cg* (thick solid 
lines), we include two charmonium vector J/^(1S) and 
4>'{2S) resonances(see Appendix A). The cusp of Re(Cg*'') 
at s = 4(mc/mb)^ ~ 0.34 as shown in Fig. |l| (thin 
line) is due to the cc-loop contribution from Y5£i(s)[see 
Eqs. (\Al\ j and (A2) in Appendix A]. In Fig. ^ one can 
also find that Rc{Cf) > ImiCf). 
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FIG. 1. The effective Wilson coefficient Cf as a function 
of s = q^/Mg. As the real part of Cf^, the thick(thin) solid 
Une represents the results with(without) LD contribution, i.e. 
Cf{Cf). The imaginary (dotted line) part of Cf is the 
result without LD contribution. 

The long-distance contribution to the exclusive B 
K decay is contained in the meson matrix elements of 
the bilinear quark currents appearing in Ti-^g^ given by 
Eq. (^. The matrix elements of the hadronic currents 
for B ^ K transition can be parametrized in terms of 
hadronic form factors as follows 

EE {K\rsj^PLb\B) = i[F+(g2)pM + F^iq')q^], (5) 

and 

EE {K\sia'"'q^PRb\B) 
1 



2{Mb + Mk) 



[q'P^ - {Mi - Mi,)q^]FT{q'), (6) 



where P — Pb + Pk and q — Pb — Pk is the four- 
momentum transfer to the lepton pair and Amf < q^ < 



{Mb - Mxf. We use the convention a^"' = (i/2)[7^, j"] 
for the antisymmetric tensor. Sometimes it is useful to 
express Eq. (g) in terms of F+{q^) and Fo(g^), which are 
related to the exchange of 1~ and 0+, respectively, and 
satisfy the following relations: 



F+(0) = Fo(0), Fo(g2) = F+(g2) + 



q 



Mi 



M^ 



F-{q'). 



(7) 



With the help of the effective Hamiltonian in Eq. (|^) 
and Eqs. (|) and (|), the transition amplitude for the 
B K£^£^ decay can be written as 



M = {K£+£-\ncs\B 
4Gf a 



V2 47r 



Gcff T 
9 '^M 



2to6 t 



(8) 



where a = e^/47r is the fine structure constant. The dif- 
ferential decay rate for the exlcusive rare B K£~^£~ 
with nonzero lepton mass(m£ 7^ 0) is given by (see Ap- 
pendix B for the detailed derivation) 



dr 



t\,t5 (-11 



3 • 297r5 



1-4 



mi 



1/2 



1 + 2-:^ \Ft+ 



s 



(9) 



where 



Fo- 



= \CfF+ - 



\c\oni-fr\Fo\ 

(s - 1 - f)^ - 4f, 



\Cio\'\F+ 



]F+\'] 



(10) 



mj/Ml, and f = Mj^/M^. 



with s — q'^/Mg, rhe ^. ^ 

We used mi, ~ Mb in derivation of Eq. (|^). Note also 
from Eqs. (H) and (|l0|) that the form factor F-{q'^){or 
Fo{q'^)) contributes only in the nonzero lepton(m£ ^ 0) 
mass hmit. Dividing Eq. (||) by the total width of the B 
meson, which is estimated to be [fZlBil 



Ftot, — 



/AfIC 



1927r 



(11) 



one can obtain the differential branching ratio dBR{B 
K£+£-)/ds = {dT{B K £+ £-)/r tot) /ds 0. 



3 and the central value of \Vcb\ = 0.0402 |3l|, we 
; 1.688 ps while r)fj? = (1.653 ± 0.028) ps. Since 



^With / 
obtain tb 

our numerical results of the branching ratios are obtained 
from using Eq. (|ll|), approximately 2% theoretical error due 
to the lifetime of B meson is understood. 
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As another interesting observable, the longitudinal lep- 
ton polarization asymmetry (LPA), is defined as 



PLi-s) = 



dVh^-i/ds - dTh=i/ds 
dTh=-i/ds + dTh^i/ds' 



(12) 



where h — +1(— 1) denotes right (left) handed I in the 
final state. From Eq. (H), one obtains for B K£^£^ 



2 1 



1/2 

4^1 ] 0CioF+ 



F+ReCl 
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(13) 

Note that our formulas for the differential decay rate in 
Eq. (^) and the LPA in Eq. (|3|) are written in terms 
of (F+ Fq Ft) instead of {F^,F^,Ft) as obtained in 
Refs. a iql. However, our formulas and those in [p|Jl0| 
are equivalent with each other once we rearrange our for- 
mulas in terms of {F+ , F_ , Ft). One nice feature of using 
-Fb in the decay rate formula is to separate the Fq contri- 
bution from the total rate as we shall show later. 



of the initial and final state light-front wave functions as 
follows 



^^oot_^A^r'^I^^oo_^ (15) 



where V = jf'PL for J'' in Eq. (||) and ia^'^q^PR for 
in Eq. (||), respectively. The measure [d^'Pq\ in Eq. (|l" 
is written in terms of light-front variables as 



dk' dk 



(16) 



where dk^/dx is the Jacobian of the variable transforma- 
tion {x, ki^ —> k — [kz, kj_) defined by 



Jo 



dk. 



dx Ax{l — x) 



1 



2 2 \ 2n 



_ml + kl ^ ml + kl 



1-x 



(17) 
(18) 



III. FORM FACTOR CALCULATION IN 
LIGHT-FRONT QUARK MODEL 

A. Analytic calculation in. — Q frame 

As shown in Eq. (|^), only two weak form factors Fj^{q^) 
and Fxiq^) are necessary for the massless(mf = 0) rare 
exclusive semileptonic b — > si^£~ process. The form fac- 
tors F+(g2) and Ft(q^) can be obtained in q+ = frame 
with the "good" component of currents, i.e. u — +, with- 
out encountering zero- mode contributions p7[ |. Thus, we 
shall perform our light-front quark model calculation in 
the q~^ = frame, where q^ = q^q^ ^'zl = < 0, and 
then analytically continue the form factors Fi(q^)(« = 
+, T) in spacelike region to the timelike > region by 
changing q± to iq± in the form factor. 

The quark momentum variables for Psiqiq) — > 
PK{q2q) transitions in the q^ —Q frame are given by 



p-t = {l-x)P+, 

Pi± = (1 - x)Pn_ + k^ 
p+ = {l-x)P+, 

P2± = (1 - X)P2J_ + k'^ 



xP+, 

Pg± = XPl± - fcj_. 



p'+ = xP+ 



2 ' 



Pg± 



= xP2\ — fc'i 



(14) 



For S ^ if 



which require that = p'^^ and Pq± — jf^^ . 
transitions, one has mi = mi,, m2 = Wg, and m^ = m,j. 
Our analysis for b s£^i~ decays will be carried out in 
this g"*" = frame and the decaying hadron (B-meson) is 
at rest, i.e. Pi_l = 0. 

The matrix elements of the currents in Eq. (|^) and 
in Eq. (^ are obtained by the convolution formula 



The spin-orbit wave function TZ'^'^'^_{x, k±) is obtained by 
the interaction-independent Melosh transformation. The 
explicit covariant form for a pseudoscalar( J = 0, = 0) 
meson is given by 



7^ 



{x, kj_) 



u{Pq,Xq)l^v{Pq,\) 

V2JmS - {ml - mlY 



, (19) 



where A's are light-front helicities. Our radial wave func- 
tion is given by the gaussian trial function for the vari- 
ational principle to the QCD-motivated effective light- 
front Hamiltonian 1201: 



{x,k±) 



1 



7r3/2/33 



1/2 



exp(-fc72/3^), (20) 



which is normalized as J d'^k\(j>{x, k±)\'^ — 1, where fc^ 
fc^ 4- fc^ and k^ is given by 



kz^{x~ i)Mo 



2 2 

2Mo 



(21) 



Then, the sum of the light-front spinors over the helic- 
ities in Eq. ( p^ are obtained as 

"a,- (P9)7^"L iP2)ux^ {p2)r^'ux, (pi)uAi (Plh^VXg (Pq) 

A's 

= Tr - m,-)75(^^2 + m2)T^'{A + ^1)7' • (22) 

Using the matrix element of the component of the 
currents(/i = -I-), and the particle on-mass shell condi- 
tion, i.e. the hght-front energy p~ = {pf^ + mf)/pf{i = 
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1, 2 and q) in Eq. (p2[), we obtain the weak form factors 
F+{q{) and Frifj) as follows 



dx 



(fk. 



A1A2 + k^-k\ 



and 



A\ 



dx 



(23) 



fc2 



d^k^ 



dx 



dx 



(/)2(a;,fc^)0i(x, /cj_) 



Mb + Mk 
(1 - x)MoM'o 



(m2 - mi) 



91 



-4i 



(24) 



where A 



xrui 



+ (1 - x)mj{i = 1,2), Mo = 
, and k', 



n,j — k± — xq±. The primed 
factors in Eqs. ( |23| ) and (p3) are the functions of fi- 
nal state momenta, e.g. k'^ = k'^{x,k'±) and M'o — 
M'o{x,k' ±). Since the weak form factors F+{q^) in 
Eq. (11) and Friqi) in Eq. (||) are defined in the 
spacelike(q^ < 0) region, we then analytically continue 
them to the timelike > region by replacing q± with 
iq± in the form factors. We describe in Appendix C our 
procedure of analytic continuation of the weak form fac- 
tors. 

Our analytic solutions will be compared with the fol- 
lowing parametric form used by many others [Q-^ 13 2^] 

^(0) 



F{q' 



1 - (71 g2 



(25) 



where the parameters cti and <T2 are determined by the 
first and second derivatives of F{q'^) at q^ = 0. 



B. Effective calculation in > frame 

Our effective calculation of weak form factors 
is performed in the purely longitudinal momentum 
frame 



IJ23] where g+ > and Pix = P2± = so 

q^q^ > is 



that the momentum transfer square q^ 
timelike. 

One can then easily obtain q^ in terms of the mo- 
mentum fraction a — P2 / Pi = 1 ^ / Pi^ — 
(1 — a){Mf — M|/a). Accordingly, the two solutions for 
a are given by 



a± 



M2 
Ml 




2M1M2 



(26) 



The -l-(— ) sign in Eq. ( p6[ ) corresponds to the daughter 
meson recoiling in the positive(negative) z-direction rel- 
ative to the parent meson. At zero recoil(q'^ = ^max) ^^'^ 
maximum recoil((j^ — 0), a± are given by 



a+{q\ 



' ) 

max/ 



,2 ^Mi 

Wmax/ ' 



M9 

a+(0) = l, a_(0)=(^ 



(27) 



The quark momentum variables in the q^ > frame are 
similar to Eq. ([l^ ) in the g+ = frame but the momen- 
tum transfer in g+ > frames flows through only lon- 
gitudinal component of quark and antiquark momenta, 
i.e. 



= {1- X)P^, = .xPi+, pi± = -Pq± = k_l 

P+ = {1- X')P+, p'+ = X'P+, P21_ = -P^q^ = , 



(28) 



where x' — x/a and Pij^ = P2±_ = has been used (see 
Fig. I). 

The a ± -independent form factors F±{q^) defined in 
g+ > frames are then obtained as follows 

P^(g^)^±(^^"-^^'^("+)-(^^"+^^'^("-\ (29) 



where j+(a±) = {K\s"f+ PLb\B)\a^/ P+ from Eq. 




rO 1-a (b) 



1-a rO 




\ 1 

FIG. 2. The covariant diagram (a) corresponds to the sum 
of the LF valence diagram (b) defined in < x < a region and 
the nonvalence diagram (c) defined in a < x < 1 region. The 
large white and black blobs at the meson-quark vertices in (b) 
and (c) represent the ordinary LF wave function and the non- 
valence wave function vertices, respectively. The small black 
box at the quark-gauge boson vertex indicates the insertion 
of the relevant Wilson operator. 

As shown in Fig. ^ the g+>0 frame requires not 
only the particle-number-conserving (valence) Fock state 
contribution in Fig. ^(b) but also the particle-number- 
nonconserving (nonvalence) Fock state contribution in 



i.e. j+{a±) =_7l,(a±) + j+„(a±) in Eq. 1^ 



Fig. 1(c); 



In our previous works I23J23], we have developed a new 
effective treatment of the non- wave- function vertex(black 
blob in Fig. ^(c)) in the nonvalence diagram arising from 
the quark-antiquark pair creation/annihilation. Since 
the detailed procedures for obtaining the effective solu- 
tion for the non-wave-function vertex have been given 
in 23|, here we briefly present the salient points of 
our effective method |22 23] and the final forms of the 



current matrix elements for both valence and nonvalence 
diagrams. 
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The essential feature of our approach is to consider 
the Ught-front wave function as the solution of light-front 
Bethe-Salpeter equation(LFBSE) given by 

[dy][(fl±_]IC{xi, hr,yj,lj±)'^{yj,lj±), (30) 



where JC is the B-S kernel which in principle includes 
all the higher Fock-state contributions, Mq — {m\ + 
klj^) / Xi + {m2 + k1^) / X2, and ^{xi, ki±) is the B-S ampli- 
tude. Both the valence(white blob) and nonvalence(black 
blob) B-S amplitudes are solutions to Eq. (|30|). For 
the normal(or valence) B-S amplitude, xi = x and 
X2 = a — X > 0, while for the nonvalence B-S ampli- 
tude, Xi — X and X2 — a — x < 0. As illustrated in 
Figs. §(b) and (c), the nonvalence B-S amplitude is an 
analytic continuation of the valence B-S amplitude. In 
the LFQM the relationship between the B-S amplitudes 
in the two regions is given by p3,E3l 



{M^ -Ml)-^'{x„h±) 

[dy] [d'^l±]IC{xi, ki±;yj,lj±)-i'{yj,lj±), 



(31) 



where (xi, ki±) represents the nonvalence B-S ampli- 
tude and again the kernel includes in principle all the 
higher Fock state contributions because all the higher 
Fock components of the bound-state are ultimately re- 
lated to the lowest Fock component with the use of the 
kernel. This is illustrated in Fig. |[ 



-X ■ 

















K 



FIG. 3. Non-wave-function vertex(black blob) linked to an 
ordinary LF wave function (white blob). 

Equations (^) and (^ are integral equations for 
which one needs nonperturbative QCD to obtain the ker- 
nel. We do not solve for the B-S amplitudes in this 
work, but a nice feature of Eq. (^) is a natural link be- 
tween nonvalence B-S amplitude and the valence one 
^ which enables an application of a light-front CQM even 
for the calculation of nonvalence contribution in Fig. ^(c). 
In (1 + 1)-QCD models ||3^,^, it is shown that expres- 
sions for the nonvalence vertex analogous to our form 
given in Eq. ( |3l| ) are obtained. With the iteration pro- 
cedure given by Eq. (|3^) in this > frame, we obtain 
the current matrix element of the nonvalence diagram in 
terms of light-front vertex function and the gauge bo- 
son vertex function. The interested reader may consult 
Refs. [HHH on this subject. 

The matrix element of the valence current, 7'^ , in 
Eq. (E9), is given by 



/ dk' dk 
"^^^-^V a^V ^'^2(a;',A:i)0i(a:,A:_L) 



+ k\ 



(32) 



Bl + kldBl + k 



where 



Bi = xmi + (1 — x)mq, B2 — x'm2 + (1 — x')mq, (33) 



and k'^ = kz{x',k±) in Eq. (|2l|). The matrix element of 
the nonvalence current, j+, in Eq. (p9|), is obtained as 



dx 



dk 



-^X^{x,k±_)4>i{x,k^) 



^ x'{l — x') J V dx 

kl + B1B2 + x{l - x){l - x'){Ml - M2) 



X J dy J d^lj 

where 

X'ix,k^) = 



^/x{l - x)Mo 

j dlz IC{x,kj_;y,lj_) 
i dy M'o(2/,4) 



9" 

1-Q 



k'j_-\-rnf k^-\-mi 
1 — x x — a 



MyJ±), (34) 



(35) 



is the light- front vertex function of a gauge boson Q and 
dy — dy/ yjy{l — y). In derivation of Eq. (^) with the 
-component of the current, we also separate the on- 
mass shell propagating part (i.e. the term proportional 
to {k\ + B1B2)) from the instantaneous part(i.e. the 
term proportional to x(l — x){l — x'){Mf — Mq)), where 
the struck quarks (mi — nib and TO2 — mg) are on-mass 
shell and the spectator quark [niq — to„) is off-mass shell. 
Note that the instantaneous contribution exists only for 
the nonvalence diagram as far as the -component of 
the current is used. As we shall show in the next numer- 
ical section, the instantaneous contribution to the weak 
form factors F±{q^) for B ^ K transition is quite sub- 
stantial near zero recoil. 

Note that Eq. ( |3l| ) was used to obtain the last term 
in Eq. (^). While the relevant operator K. is in general 
dependent on all internal momenta (x, k±;y,l±), the in- 
tegral of /C over y and l± in Eq. (^4|) depends only on x 
and k±, which we define 



GBKix,k±) = dy d^li 



I dlz IC{x,k±;y,l±) 
dy M'o(y,/l) 



4>2{y,l±)- 
(36) 



■^While one can in principle also consider the B-S amplitude 
for we note that such extension does not alter our results 
within our approximation in this work because both hadron 
and gauge boson should share the same kernel. 
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In this work, we approximate GBK{x,k±) as a constant 
which has been tested in our previous works ||2^ , p3|] and 
proved to be a good approximation. As we shall show in 
the next section, the reliability of this approximation can 
be checked by examining the frame-independence of our 
numerical results. 



IV. NUMERICAL RESULTS 

In our numerical calculation for the process of i? — > 
K£^£~ transition, we use the linear potential parame- 
ters presented in Ref. . Our predictions of the decay 
constants for K and B were reported [^,^ as //f =161.4 
MeV(Exp.= 159.8±1.4) |§ and fs = 171.4 MeV [H, 
respectively.^ Our model parameters and decay con- 
stants are summarized in Table | and compared with ex- 
perimental data 1^ as well as lattice results |3^. Note 
that in the numerical calculations we take mf, — 5.2 GeV 
in all formulas except in the Wilson coefhcient Cg^ , where 
?7ifc = 4.8 GeV has been commonly used. 



{Mb — MkY ■ We also include the results obtained from 



2.0 
1.5 
1.0 
o 0.5 



|i 0.0 




-1.0 
-1.5 
-2.0 



Fjq ) in q =0 frame 
F^(q^) in q'=0 frame 
F{q') fromF(0)/[1-o,q'+Ojq'] 
F,(qVromF,(0)/[1-o,q%a,q'] 



5 10 

q'[GeV'] 



15 



20 



FIG. 4. Analytic solutions of (g'^) (thick solid line) and 
FT(q'^)(thick dashed line) compared with the results(thin 
lines) obtained from the parametric formula given by Eq. ( |2^ ) 
for B ^ K transition. 

In Fig. ^ we show our analytic (g"*" = frame) solu- 
tions for the weak form factors F+(g^) (thick solid line) 
and Ft (q^) (thick dashed line) for -5 GeV^ < < 



The difference of decay constants between this work and 
Refs. po[pl| is only due to the definition, i.e. we use the def- 
inition (fi\q2Yl5qi\P) = ifpP'^ in this work so that /^''p- = 
130.7±0.1 MeV while we used (0|g27^75gi|P) = iV^fpP" in 
Refs. EoiP 



the parametric formula given by Eq. ( |25D where the thin 
solid (dashed) line represents F+( Ft)- Our analytic solu- 
tions given by Eqs. (E3|) and (|24| ) are well approximated 
by Eq. ( p5| ) up to (T ^ 15 GeV^ but show some de- 
viations near zero recoil point. We summarize in Ta- 
ble p our numerical results for the weak form factors 
F_|_(g^) and FtIq'^) at = and the parameters ai 
defined in Eq. (|25|) and compare with other theoreti- 
cal results ||7|,P, |13| , p9| . As one can see from Table ||, 
our results for the F^{q^) and FT{q^) in g-^ limit 
are quite comparable with other theoretical results. As 
other theoretical schemes predicted, our results also show 
F+(0)(= 0.348) ~ -Ft(0)(= -0.324). 



2.5 
2.0 
1.5 
1.0 

0.5 
0.0 
-0.5 
-1.0 



Anaiytic soi. from q*=0 frame 
Effective sol.(=val.+nv.) from q*>0 frame 
Valence part from q*>0 frame 
Instantaneous part from q*>0 frame 




10 15 
q'[GeV'] 



20 



FIG. 5. Effective solution of F+{q^){t\an solid line) for 
B ^ K transition. The line code is in the figure. 

For the analysis of heavy r decay process, the weak 
form factor F_((7^)(or equivalently Fo(g^)) is necessary 
for the calculations of the decay rate and the LPA and 
we obtain it using our effective method (22 j2^ in (7+ > 
frame as described in Sec. III(B). In FigT^, we show our 
effective((7"'" > frame) solution of F+{g^) (thin solid 
line) with a constant Gbk = 3.9 fixed by the normal- 
ization of F+{q^) in the (7+ = frame (thick solid line) 
aX q-^ — limit. As one can see in Fig. |[ our effec- 
tive solution of F+((7^)(thin solid line) is very close to 
the analytic one (thick solid line) for the entire kinematic 
region. It justifies the reliability of our constant approx- 
imation Gbk of the kernel K,. For comparison, we also 
show the valence(dottcd line) and the instantaneous (dot- 
dashed line) contributions to Fj^{q^) in the g"*" > frame. 
Although the valence contribution dominates over the 
nonvalence one for q^ < 10 GeV^, the nonvalence (es- 
pecially the instantaneous) contribution is not negligible 
for g2 > 10 GeV2. 
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1.0 



Effective sol. from q >0 frame 
Valence part from q*>0 frame 
Instantaneous part from q'>0 frame 



0.4 



presented in Fig. ^ agrees very well with the light 
cone QCD sum rule (LCSR) result for F_[q^) by Ahev 
et al. HKSee their Fig. 1(b)). Similarly, our effective so- 
lution for ^0(9^) is in a close agreement with the LCSR 
results given by Ball jl^ and Ali et al. jl^. Our effec- 
tive solution of Fo(9^) as well as the analytic solutions 
of F+{q^) and Friq^) shown in Fig. ^ will be used for 
the calculations of the branching ratios and the longi- 
tudinal lepton polarization asymmetries. We shall also 
discuss how we take the effect of the vector meson dom- 
inance (VMD) into account at the end of this section. 



0.2 - 



0.0 



10 15 
q'[GeV'] 



20 
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resonance + nr 
nonresonance(nr) 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.f 



q/Mg 



FIG. 6. Effective solutions (solid line) of Fo(g^) and F^{q^) 
compared with tlie valence contributions(dotted line) for 
B ^ K transition. 

Using the same constant operator Gbk = 3.9, we 
are now able to calculate the scalar form factors F^iff) 
and F^(q^) in (7+ > frames and the results are 
shown in Fig. ^(solid line). As in the case of F+{q'^) in 
Fig. ^, we also include the valence contributions (dotted 
line) to both Fo{q^) and F-.{q^) and the instantaneous 
contribution(dot-dashed line) to ^0(9^)- It is very in- 
teresting to note especially from F_(<j^) that the non- 
valence contribution, i.e. the difference between solid 
and dotted lines, is very substantial even at the max- 
imum recoil point((7^ = 0) and is growing as q^ in- 
creases. As a reference, our numerical results for F^ ob- 
tained from our effective (valence) solution at maximum- 
and zero-recoil limits are F_(0) = — 0.14(— 0.34) and 
-^-(^max) = — 0.9(— 2.23), respectively. Our result for 
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FIG. 7. The branching ratios for B Ki+T = e, /i) (a) 
and B Kt^t~ (h) transitions. The thick(thin) solid line 
represents the result with(without) LD contribution to C9* 



in Eq. 



The dotted line in (b) represents the Fo(q ) con- 



tribution to the total branching ratio of r decay. 
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We now show our results for the differential branch- 
ing ratios for B K£^£~{£ = e,/i) in Fig. 0(a) and 
B Kt^t^ in Fig. 0(b), respectively. The thick(thin) 
solid line represents the result with(without) the LD 
contribution(YL£)(s)) to Cg^ given by Eq. (|). In plot- 
ting Figs. 0(a) and (b), we set nii = and nir — l-lll 
GeV, respectively. As one can see the pole contributions 
clearly overwhelm the branching ratio near J/ip(lS) and 
tp'i^S) peaks, however, suitable £~^£~ invariant mass cuts 
can separate the LD contribution from the SD one away 
from these peaks. This divides the spectrum into two dis- 



tinct regions 



(i) low-dilepton mass, imj 



< 



< 



S, and (ii) high-dilepton mass. Mi, + S < q < 



9max' where 6 is to be matched to an experimental cut. 
The branching ratios with[without] the pole(i.e. LD) con- 
tributions for B Kl'^£~ are presented in Table H for 
low(second column), high(third column), and totaI(4th 
column) dilepton mass regions of . Although the con- 
tribution of scalar form factor F^^q^) to massless lepton 
decay is negligible(zero for mi = 0), its contribution to r- 
decay as shown in Fig. 0(b) (dotted line) is very substan- 
tial, e.g. ^ 75% contribution to the total(nonresonant) 
decay rate in our model calculation. Thus, the reli- 
able calculation of FQ{q^) is absolutely necessary and our 
effective method of calculating the nonvalence diagram 
seems very useful. 

It is worthwhile to compare our results for the branch- 
ing ratios with other light-front quark models |p|,p^. 
While the authors in Ref. used the simple parametric 
formula, Eq. (p5|), to obtain F+ and Ft and the heavy 
quark symmetry(IIQS) to extract F^ , the authors in 
Ref. used the dispersion representation through the 
(Gaussian) wave functions of the initial and final mesons 
and then analytically continue the form factors from the 
spacelike region to the timelike region. The common as- 
pect in these models is to have the same form factors 
and i^Tj which are free from the zero-mode contribution, 
not in the timelike region but in the spacelike region as 
far as the same model parameters are used. Indeed our 
method of analytic continuation of the form factors F^ 
and Ft is equivalent to that of Ref. |0 . However, the dif- 
ference is in the calculation of F_ , which is not immune to 
the zero-mode contribution. The zero-mode contribution 
must be properly taken into account for the calculation 
of F^ . Thus, it is not quite surprising to note that al- 
though our branching ratio(see Fig. 0(a)) for the massless 
lepton [£ = Cjji) decay is not much different from the re- 
sults in Ref. [g[(see their Fig. 1(a)) and Ref. jlO|(see their 
Fig. 3(a)), our branching ratio(see Fig. 0(b)) for the r de- 
cay is quite different from the results in Ref. Q (see their 
Fig. 1(b)) and Ref. 0(see their Fig. 3(c)). 

Our numerical results for the non-resonant branch- 
ing ratios (assuming \Vtb\ — 1) are 4.96 x 10~''| Vts/KfeP 
for B Kl+£- {£ = e,fi) and 1.27 x 10-'^\Vts/Vcb\'^ 
for B Kt+t^, respectively. While the CLEO Col- 
laboration |l[] reported the branching ratio Bt(B 
Ke+e-) < 1.7 x 10~^ the Belle Collaboration(K. Abe 



et al.) reported Br(B Ke+e') < 1.2 x 10"^ and 
Bt{B Kii+^r) ^ (0.99t°j^t°/i^) X 10-^ respec- 
tively. Our non-resonant results for the branching ratios 
oi B ^ K£^£^ is summarized in Table IV and compared 
with experimental data as well as other theoretical pre- 
dictions within the SM. 

The exclusive B — > Kt^t~ has been computed via the 
heavy meson chiral perturbation theory by Du et al. [p^ , 
where the branching ratio of the exclusive decay was 
found to be about 50—60% of the inclusive one. Although 
calculations of exclusive decay rates are inherently model 
dependent, chiral perturbation theory is known to be 
reliable at energy scales smaller than the typical scale 
of chiral symmetry breaking, Acsb — 47r/^/-\/2. In 
B Kt^t~, the maximum energy of the fC-meson 



in the B rest frame is {Mg 



Ml 



'iml)/2ME 



1.5 



GeV, which places most of the available phase space 
around the scale Acsb 24|. From the above argu- 
ment and our exclusive t branching fraction, we can es- 
timate the branching ratio of inclusive B — > XsT^t^ 
as (2.12 - 2.54) x lQ~'^\Vts/Vcb? which is quite com- 
parable to the prediction given by Hewett |Q where 
BR(B XsT+T-) = 2.5 x 10-"^ was obtained. 

In Figs. ||(a) and (b), we show the longitudinal lepton 
polarization asymmetries for B — > Kfi^fj,^ and B — > 
Kt^t~ as a function of s, respectively, and with (thick 
solid line) and without (thin solid line) LD contributions. 
For the B Kfi^n^ case, we use the physical muon 
mass, r7i^=105 MeV. In both figures, the longitudinal 
lepton polarization asymmetries become zero at the end 
point regions of s. Our numerical values of Pl without 
LD contributions and away from the end point regions 
are -0.97 < Pl < -0.98 in 0.3 < s < 0.6 region for 
B Kfi+fi" and -0.15 < Pl < -0.18 in 0.5 < s < 0.7 
region for B Kt'^t~ , respectively. In fact, the Pl 
for the muon decay is insensitive to the form factors, e.g. 
our Pl — — 0.98(away from the end points region) is well 
approximated by |11[ 



Pl 



CinReCS*^ 



\Ci 



I. 



(37) 



in the limit of Cy — > from Eq. (|l3|) . It also shows that 
the Pl for the /i dilepton channel is insensitive to the 
little variation of C7 as expected. On the other hand, 
the LPA for the t dilepton channel is sensitive to the 
form factors. In other words, as in the case of branching 
ratios, although our result of the LPA for the muon de- 
cay is not much different from the results in Ref. [^(see 
their Fig. 2(a)) and Ref [0(see their Fig. 5(a)), the 
result for the tau decay is quite different from the re- 
sults in Ref. §(see their Fig. 2(b)) and Ref 0(see their 
Fig. 5(c)). 
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FIG. 8. The longitudinal lepton polarization asymmetries 
Pl(s) for B Kl+r{a.) and B KT+T-{h) transitions. 
The same line code is used as in Fig. ^. 

Comparing our results for the weak form factors with 
other phenomenological models, one may find that there 
is in general a good agreement for small and intermedi- 
ate region. Nevertheless, there are some differences for 
large region where vector mesons are expected to dom- 
inate(VMD) especially for For example, both 

results of the LCSR in [|l3|,|9) and our LFQM analyses 
show that the direct solution for F^{q'^) is well approx- 
imated by Eq. up to q^ < 15 GeV^. However, the 
large momentum behavior of F_|_(g^)(as well as Fxiq'^)) 
is somewhat different since our model does not include 
the VMD effect. 

Following the same method used in recent LCSR anal- 
ysis we use the VMD formula(i.e. i?*-pole with 
Mb- = 5.325 GeV) given by 



1 - q^/Ml 



(38) 



at large region and match the parametric formula 
in Eq. (p5|) by the following constraint |59| 



F^^D(<Z^)=F+(g^)inEq. (25), 



dq2 



VMD 



(9') = ^i^+(9')inEq. (|2 



(39) 



to make both parametrizations smooth connection at a 
transition point q^ =90' where c is fixed at q^ = 9q ™ 
Eq. (H). We should note that the F+{q^) in Eq. (^) is 
almost equivalent to our LFQM prediction F]f^^{q^) up 
to q^ ^ 15 GeV^ and the transition point is expected 
to be at q^ ~ 15 GeV2(see also Ref. [M) in order to make 
interpolation between F^^'^^[q'^ and F^^^{q^ > 

Qq) more sense. ^ In our case for B K transition, we 
obtain (c, (7§)=(0.388, 14.38 GeV^) for Ff^{q^). For the 
tensor form factor, we get (c, gg)=(-0.358, 14.23 GeV^) 
for Friq'^). 

It is necessary to discuss the exclusive B irivi pro- 
cess in that the constant c has a direct physical implica- 
tion for B — > Tr£i^i process, i.e. it is related to the physical 
couphngs as gEfLpl 



2Mb' 



(40) 



where Jb* is the decay constant of the B* meson defined 
by {0\b'y''u\B*) = Mb'/b'^'" and 35.5^ is the (axial- 
current) coupling defined by {B°{P)TT+{q)\B*+{P + 
q)) = gBB-'Tr{q ■ e) and can be extracted from soft pion 
q^ limit in the heavy meson chiral perturbation 
theory 44|. In the limit where the heavy quark 
mass mQ[Q = c, h) goes to infinity there are flavor- 
independent relations between coupling constants 



g 



2Mt 



-go'D-iT 



2Mf 



-gB'B-jT, 



(41) 



where /^r = 131 MeV and the coupling constant g appears 
in the interaction Lagrangian of the effective meson field 
theory 0,11,0. 

In our numerical calculation of c for the exclusive 
B Tieve process, we obtain (c, g§) = (0.312,15.12 GeV^) 
from E q (|3g| ) and (cri,CT2) = (4.75 x 10-2,5.50 x 10"") 
in E g. ( |2^ ) , which was obtained in our previous anal- 
ysis [Q. Since we also obtained the B* meson decay 
constant as /s* = 185.8 MeV we can now ex- 

tract the coupling constant of the B* to STr-pair and 



''As discussed in [^, a naive extrapolation of the VMD 
formula in Eq. ( |38| ) to the point = is not consistent with 
the monopole formula F+{q^) = F+(0)/(l - q^/Al) used in 
many theoretical ansatz since the relevant parameters are in 
general different, i.e. -F+(0) 7^ c and Ai 7^ Mb* ■ 
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the result is gB*B-K+ = 17.88 and g=0.23 while the re- 
cent fit [Em to the experimental data gives two possible 



solutions, g — 0.27 



+0.04H 
-0.02- 



0.05 
0.02 



or O - I) Yh W.03-F0.2 ^ 
m y — ^- '"-0.03-0.1- 

acknowledge the remark in that for the B -kIvi 
form factors with E'^ < 2rn^, analytic bounds combined 
with chiral perturbation theory give g/s % 50 MeV . 
That means while the solution g — 0.27 gives /s < 190 
MeV, g = 0.76 gives /b ^ 66 MeV, which is roughly a 
factor of three smaller than lattice QCD result |^,i.e. 
/Lat- = 200 ± 30 MeV. Note that our LFQM prediction 
is given by j^'^^ — \1\A MeV. As a reference, other 
theoretical calculations for g are 0.2 — 0.4 for the QCD 
sum rules, 1/3 — 0.6 for the quark models^ and 0.42(4)(8) 
for the lattice calculation(see Ref. [Q for the survey of 
g values obtained from different models). 



VMD CoiTCction to LFQM: 
F^(0)/ri- <J,q' + O,q^l(q'=[0.q„-1) 
- c/[l-q7m\.Kq'=rq„'.q'„,„l) 




FIG. 9. VMD corrections to the LFQM predc- 

tions for ^(g2)(solid line) and ^^''(g^) (dashed line), i.e. 



In Fig. ^, we show the VMD corrections to both 
""(q^) (solid line) and F_P^(g^) (dashed line), i.e. 
F+(g2) = F_^:^Q%2 < ^2)^^MD(^2 > ^2)^ Comparing 
Fig. H[Fig. 3 in [pl|] and Fig. 0, we find the enhancement 
of Ff ^(g2)[i?S'^(g2-)] at q2 = ql^^ by around 40[70]%. 

Our result for Ff'^{q^) including the VMD correction 



pBK 



^Using similar LFQM to ours, Jaus obtained g = 0.56 
from the direct calculation of the hadronic matrix element in 
the soft pion limit and argued that the calculated p — -k — n 
and K* —K — n coupling constants within the same model are 
in fair agreement with data. The reason for the discrepancy of 
g value is not yet understood. However, the computed decay 
constants /s and fs* are in good agreement between Ref. 
[Eol and ours. 



are quite coraparable with that obtained from QCD sum 
rules in Ref. |39(| where the authors used the same method 
to enhance F^(g^). Our result for -F£f"(V) in Fig. ^ 
is also comparable with those of Refs. |l^,|l^. However, 
the branching ratio for B K£^£^{£ — e^fi) increases 
less than 2% by including the VMD effect. It is not sur- 
prising to note that the large enhancement of the weak 
form factors near the zero-recoil(g^ = q^ax) region does 
not affect the differential decay rate very much, since the 
phase space of the large q^ region is highly suppressed in 
Eq. (P). 



V. SUMMARY AND CONCLUSION 

In this work, we investigated the rare exclusive 
semilpetonic B K£^£^ {£ — e, jj, and r) decays within 
the SM, using our LFQM which has been tested ex- 
tensively in spacelike processes [^,^ as well as in the 
timelike exclusive semileptonice decays of pseudoscalar 
mesons [|l|,||. The form factors F+{q'^) and Friq"^) 
are obtained in the q^ — frame (g^ < 0) and then 
analytically continued to the timelike region by chang- 
ing q± to iq±^ in the form factors. The form factor 
F-{q^) is obtained from our effective treatment of the 
nonvalence contribution in addition to the valence one 
in 9+ > frames {q^ > 0) based on the B-S formal- 
ism. The covariance (i.e. frame-independence) of our 
model has been checked by comparison of F+{q'^) ob- 
tained from both q~^ — and g"*" > frames. Our nu- 
merical results for the form factors are comparable with 
other theoretical calculations as shown in Table Us- 
ing the solutions of F+ and Ft obtained from g+ = 
frame and F_ obtained from g+ > frame, we calcu- 
late the branching ratios and the longitudinal lepton po- 
larization asymmetries for B K£^£^ including both 
short- and long-distance contributions from QCD Wilson 
coefficients. Our numerical results for the non-resonant 
branching ratios are in the order of 10~^, which are con- 
sistent with many other theoretical predictions as shown 
in Table IV. Of particular interest, we were able to esti- 
mate the inclusive branching ratio for B X^t^t' as 



BR(B ^ X,T+T-) ~ (2.12-2.54) x lQ-'^\Vts/Vcb? with 
the help of chiral perturbation theory |l^ . For the LPA 
as a parity-violating observable, we find that the LPA for 
the T channel is sensitive to the form factors while the 
LPA for the /i channel is insensitve to the model for the 
hadronic form factors. Thus, the experimental data of 
the LPA for r decay would provide a useful guidance for 
the model building of hadrons and make a definitive test 
on existing models. 
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APPENDIX A: FUNCTIONS Ysb{s), Yi,o{s), 
AND uj{s) in Eq. (0) 



The function 1sd(s) in Eq. is given by 
Fsd(s) =/i(m„s)(3Ci+C2+CW) 



2 
1 

2 

v:sVub 



Ce) 



(3Ci 



3C4 



C2)[hi0.,s)-h{m„s)], (Al) 



where C(°) = 3C3 + C4 + SCg + Cg. The function 
hlrhq—mq/mb^s) in Eq. (Al) arises from the one loop 
contributions of the four quark operators Oi — Og and 
h{mc, s), h{l, s), and /i(0, s) represent c quark, b quark, 
and u, d, s quark loop contributions, respectively. The 
explicit form of h{rhq, s) is given by 



^Inm, 



h[mq,s) = "gli^l ~ 

^i2 + yq)^\l-yq 



27 ^ 9^" 



9 



In- 



x\e{l-yq) 

~Q{yq ~ l)2arctan 



'I 


-Vq 


l\ 


-Vq 







\jyq 



(A2) 



adopt K=2.3 to reproduce the rate of decay chain 

B XsJ/^ xj+e-. 

In the SD contribution of 6 — s- s£^£^ , the u-quark loop 
contribution is neglected due to the smallness of the con- 
tribution Kb/ V;*ytfc ~ 0(A2) (A ~ 0.22 is Wolfenstein 
parameter) compared with V*gVcb — — Vj*Vt6. The term 
{V*sVub/VflVtb)C^°^ in LD contribution is also neglected 
ioi b ^ si+i-. 

The function a;(s) in Eq. represents the 0{as) cor- 
rection from the one-gluon exchange in the matrix ele- 
ment of O9 pi: 



uj(s) 



9 3 

5 + AS 
"3(1 + 2J) 

5 + 9s- 6s^ 



4 2 
Li2(s) - -lnsln(l - s) 



3(l-s)2(l + 2s) 



6(l^.s)(l + 2s)' 
where Li2{x) = — dt ln(l — xt)/t. 



(A5) 



APPENDIX B: DERIVATION OF THE 
DECAY RATE FOR B !<£+£- 

In this appendix, we show the derivation of the decay 
rate for B — > K£^£^. For simplicity, we shall omit the 
factor V(* Vtf, in the following derivation. 

The transition amplitude for B — > K£^£^ is given by 



M = {K£+£-\n\B) 



AGf a 
^/2 An 



Cfj, 



2mh rj. 



(Bl) 



where yq — Arhq/s and 

/i(0, s) — ln[ — -\ Ins H — iir. 



27 9 V M , 
The function 1ld(s) in Eq. (0) is given by 



. , 3k 



(A3) 



VtlVtb 



E 

Vi = .J/^,ip', 



TTr{V, £+l-)Mv, 
- sml - iMvTv, ' 



(A4) 



where r{Vi £~^£~), Ty. and My. are the leptonic 
decay rate, width and mass of the ith 1 cc reso- 
nance, respectively. In our numerical calculations, we use 



T{J/iP £+£-) = 5.26 X 10-6 QeV, M, 



3.1 GeV, 



Tji^ = 87 X 10-*^ GeV for J/V'(1S') and T{iP' £+£-) = 
2.12 X 10-6 GeViVV = 3.69 GeV, F^^ = 277 x 10-^ 
GeV for ^jj'{2S) |31[. The fudge factor K is introduced in 
Eq. (A4) to account for inadequacies of the naive fac- 
torization framework (see 1321 for more details.) We 



For all possible spin configurations, we make the replace- 
ment 



\M\' ^ \M\-^ 



1 



(25s + 1)(25k + 1) 



all spin states 



(B2) 



where SsiSx) is the spin of B{K) meson and we sum 
over the spins of the lepton pair. After summing over all 
spin states for the lepton pair, we obtain 



2^2" 



{2{P ■ Vi){P ■ VI) 



2J1 



T+ 



s 



(B3) 



where -F't+ is given by Eq. (nO) and 



.Fo+ = ICiol' ( b'P' - (P • qf\\F^\^ + (P • qf\F, '2 



(B4) 
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Here, we use mt ~ Mb in the deriv atio n of Eq. 

In the _B— meson rest frame, Eq. (B3) can be rewritten 

as 



[\pKf - {Ei - Eef]FT^ 



(B5) 



where \Pk\'^ = Ml(j)/4. 

The differential decay rate for B Kt^t'^ is given by 



dT 



d'pK 



2Mb \ (2T:f2EKj \ {2T:f2Ei 
■x{2^)H\Pb~Pk-Pi~Pi) 



(27r)32i?,- 



(B6) 



After doing the P^ integration, one obtains 



dT ^ 



MbGp 2 

-IT-. — « 

647r5 



[iPxf - {2Ei + Ek - MB f]FT^ 



dEKdEi 



(B7) 



The lepton energy i?^ in Eq. (Bl) satisfies the following 
upper (i?^) and lower (iJ^) bounds 



P± {MB-EK)±\PKWl-A{m,/s) 



Finally, the integration of Eq. (B7) over Eg with dEx — 
{MB/2)ds gives Eq. (|). 

APPENDIX C: ANALYTIC FORM OF THE 
WEAK FORM FACTORS IN TIMELIKE 
REGION 

In this appendix, we show the generic form of our analytic 
solutions for the weak form factors F+((7^)[Eq. (|2^)] and 
JT((7^)[Eq. (24)] in timelike region. 

In our numerical analysis, we use change of variables 

as 



(Cl) 



Since the form factors in Eqs. ( p3[ ) and ( |24D involve the 
terms proportional to {t±^ ■ which are related to 

the imaginary parts of the form factors by changing q±^ to 
«g_L, we separate the terms with even powers of (f^) 

from those with (£x • 'fj_)°'^"^ in the form factors. One 
useful identity in this separation procedure is 



V2^a + b{p^ ■q±) = \/a+^/^P{p± ■ qx_f 



(C2) 



' q"^ cos t 



By changing pj_ ■ qi_ ip± ■ q± = i\e±\ 
iSi where q^ > 0, we separate the 'Real'-parts from 
'Imaginary'-parts in Eqs. (ES) and (pih as follows 



iR{ii,q') + i6iei{ii,q'), 



from the exponent of and 

JB{i'^,q^) + i5iJi{i'^,q^), 



dx 



dx 



(C3) 



(C4) 



from the Jacobi factor. The separations of Eqs. ( |C3| ) 
and (C4) are common for both and FT{q^)- The 

main difference between the two form factors comes from 
different vertex structure and we denote generically as 



AiA 



X's 



MR(li,q 



i5iMi{Ii,q^) 



(C5) 



Combining Eqs. (C3-C5), we separate the 'Real' and 
'Imaginary' parts of the weak form factors: 



F{q' 



1 



dx I exp(— 



(7r/3i/32)3/2 7o 
[JrMr - 6'iJiMi][cos{5iIi) ~isin{5ili)] 

+5i [JrMi + JiMr] \sin{5iii) + « cos((5,l/)] , 
FflXg') + iF,,„(g2). (C6) 



We do not list here the detailed functional forms of other 
terms. However, since only the term 5i is of odd power 
in and (f^, one can easily check the imaginary term 
of the form factor Fjjn{q'^) vanishes after £± integration 
due to the fact that / d^ejL£°/'^_exp{-e''^'"^) = 0. In fact, 
we also found that the term Silj is small enough to make 
cos((5;£/) ~ 1 and sin((5;£/) ~ Silj with very high accu- 
racy. 
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TABLE I. Model Parameters (m,,/3) and the decay con- 
stants defined by (0|g27''75Qi|-P) = ifpP^ for tt, K and B 
mesons used in our analysis. We also compare our decay con- 
stants with the data jsij and the lattice result |37|. 



Meson{qQ) 


mQ[GeV] 


/3,Q[GeV] 


/ [MeV] 


ycxp. 




0.22 


0.3659 


130 


131 


K 


0.45 


0.3886 


161.4 


159.8±1.4 


B 


5.2 


0.5266 


171.4 


200± 30 Q 


TABLE II. Results for form factors F{Q) and parameters 
CTi defined in Eq. (|25|). 


Model 


F+(0) (Tl (T2 


Ft (0) (71 (72 


This work 


0.348 4.60E-2 5.00E-4 


-0.324 4.52E-2 4.66E-4 


QM 


l\ 


0.30 6.07E-2 1.08E-3 


-0.30 6.01E-2 1.09E-3 


QMJ 


9| 


0.36 4.8E-2 6.3E-4 


-0.346 4.9E-2 6.4E-4 


SR [1 




0.341 5.06E-2 5.22E-4 




SR 1 


1 


0.35 4.91E-2 4.50E-4 


-0.39 4.91E-2 4.76E-4 



TABLE HI. Branching ratio(in units of \Vts/Vcb?) 
with [without] the pole contributions for B Kt^l" for low, 
high, and total dilepton mass region. 



Mode 


1 < g2 < 8 


16.5 < < 22.9 


4m| < g^ < 22.9 [GeV^ 


(e,^i) 


2.59 X 10"'' 


3.34 X 10"** 






[2.25 X lO"'^] 


[3.70 X 10"**] 


[4.96 X lO"''] 


r 




7.20 X 10"** 
[7.47 X 10"**] 


[1.27 X 10"'^] 
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TABLE IV. Non-resonant branching ratio(in units of 
10"'^ X \Vts/Vcb\^) for B Kl+r transition compared with 
other theoretical model predictions within the SM as well as 
the experimental data taken from the Belle Collaboration (by 
K. Abe et al.) il. 



Mode This work [L0| 1 



1161 



Exp. 



e 
r 



4.96 
4.96 
1.27 



4.4 
4.4 
1.0 



3.2 ±0.8 
3.2 ±0.8 
1.77 ±0.40 



5.7 
5.7 
1.3 



< 1.2 X 10"'' 

(o.99«:it°:f5) X 10- 
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